General relativistic analysis of peculiar velocities 



George F. R. Ellis 1 Henk van Elst 2 | and Roy Maartens 3 | 



1 Cosmology Group, Department of Mathematics and Applied Mathematics, University of Cape 
Town, Rondebosch 7701, South Africa 

2 Astronomy Unit, Queen Mary, University of London, Mile End Road, London El 4NS, United 
Kingdom 

3 Relativity and Cosmology Group, School of Computer Science and Mathematics, Portsmouth 
f - *) University, Portsmouth POl 2EG, Britain 

o ' 

Abstract. We give a careful general relativistic and (l+3)-covariant analysis of cosmological 
f) |. peculiar velocities induced by matter density perturbations in the presence of a cosmological 

, constant. In our quasi-Newtonian approach, constraint equations arise to maintain zero shear of 

i the non-comoving fundamental worldlines which define a Newtonian-like frame, and these lead 

to the (l+3)-covariant dynamical equations, including a generalized Poisson-type equation. We 
investigate the relation between peculiar velocity and peculiar acceleration, finding the conditions 
under which they are aligned. In this case we find (l+3)-covariant relativistic generalizations of 
well-known Newtonian results. 
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1. INTRODUCTION 

Density and velocity perturbations of cold matter models of the observable part of the expanding 
(3JT), Universe are central to the analysis of structure formation in cosmology (see e.g. the recent review 

by Dekel ||). On scales well below the Hubble radius, the Newtonian theory of gravitation is a 
good approximation, and Peebles' approach [^lj is widely used (see e.g. [[l5| |[ [l9[ ^0|). However, 
■ observations and simulations are probing scales which are a significant fraction of the Hubble radius, 

thus requiring a relativistic treatment of the gravitational dynamics in a cosmological context. 
Furthermore, observational evidence for a positive cosmological constant (see e.g. Perlmutter 
et al ]2^]) implies that this constant does affect the evolution of matter density and velocity 
perturbations on sufficiently large scales. 

Relativistic gauge-invariant perturbations of Friedmann-Lemaitre models are usually studied 
via a metric-based approach . Bardeen's quasi-Newtonian (or longitudinal) gauge is well adapted 
to study the evolution of matter perturbations 0, [l^, 0J. Such an approach has recently been 
used to find relativistic corrections to the Newtonian relation between the peculiar velocities 
of matter structures and the matter density contrast []l7|, |2^] . Here we follow an alternative 
relativistic approach, based on the (l+3)-covariant analogue of the quasi-Newtonian gauge, which 
was developed in E3|. The quasi- Newtonian gauge is based on time slices with zero-shear normals. 
The (l+3)-covariant analogue uses a 4-velocity field u a (or a threading of spacetime) which is 
irrotational and shcarfree, and so effectively mimics the gauge fixing conditions of the quasi- 
Newtonian gauge. Analysis of the evolution and constraint equations shows that consistency 
conditions arise from enforcing vanishing shear and vorticity []To|| . These consistency conditions 
lead to the perturbation equations. 
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Vanishing vorticity reflects the absence of vector modes, while vanishing vorticity and shear 
together imply vanishing gravito-magnetic Weyl curvature, so that tensor modes are also absent. 
The remaining, scalar, modes describe peculiar velocity and matter density perturbations. Peculiar 
velocity in the matter distribution is directly generated by inhomogeneity in the expansion rate of 
u a , which also accounts for time-dependence of the peculiar gravitational potential. The peculiar 
velocity and acceleration satisfy a coupled pair of evolution equations. The peculiar gravitational 
potential obeys a generalized Poisson-type equation, which on small scales reduces to the spatial 
gradient of the usual Poisson equation. 

The peculiar velocity is not in general aligned with the peculiar acceleration, but its orthogonal 
component will remain zero if it vanishes initially. In general, the magnitude of the orthogonal 
component decays with expansion, but in order to keep this component strictly zero, we arrive at 
an equation determining the ratio of peculiar velocity to peculiar acceleration. This leads to (1+3)- 
covariant relativistic generalizations of well-known Newtonian relations, including a cosmological 
constant. 



2. (l+3)-COVARIANT QUASI-NEWTONIAN ANALYSIS 

Peculiar velocities must be defined relative to a preferred frame with 4- velocity field u a . This 
preferred frame is necessarily non-comoving relative to the matter, since, by definition, peculiar 
velocities vanish in the comoving (Lagrangian) frame. (Much of the standard literature on structure 
formation uses the term "comoving" in a rather misleading fashion, i.e. as comoving with the 
fictitious background 4-velocity rather than with the actual matter 4-velocity.) A (l+3)-covariant 
(and local) approach seeks to define the preferred frame in a physical way. (For a recent review 
on (l+3)-covariant methods see e.g. ||; for a detailed discussion of the relation between the 
(l+3)-covariant and metric-based approaches see |B.) Motivated by the Newtonian analysis of 
Peebles jnj, we choose u a to be irrotational (u> a — 0) and shearfree (a a b = 0), defining a quasi- 
Newtonian frame [jioft corresponding to Bardeen's quasi-Newtonian gauge 0. Furthermore, we 
require that u a is non-relativistic relative to observers comoving with the matter, whose 4-velocity 
field u a is given by 

u a = u a + v a , v a u a = , v a v a < 1 . (1) 

Here v a is the non-relativistic peculiar velocity, which vanishes in the background (and hence is 
gauge-invariant) . 

Thus we have two physical frames, defined by the (pressure-free) matter and by the quasi- 
Newtonian observers. In the comoving (Lagrangian) frame ii a , the pressure, energy current density, 
anisotropic pressure and 4-acceleration all vanish: 

P = <f = tab = o , A a = o . (2) 

In the quasi-Newtonian (Eulerian) frame u a , 

uj a = <r ab = 0. (3) 

The Galilean velocity boost in ([!]) preserves the pressure and energy density to linear order in v a , 
i.e., p = p, p, = fx, but introduces a non-zero energy current density (momentum density) in the 
Eulerian frame given by 

q a = pv a , (4) 

where p is the background matter density. Also to linear order in v a we have ir a b = 0. (These 
linear results are obtained from the exact boost transformation laws (B7) - (BIO) given in [fl6|| .) 
In addition, it transforms the kinematic quantities to linear order in v a as follows (cf. the exact 
boost transformation laws (B3) - (B6) in [p^| ): 

9 =0 + D a v a , (5) 
A a = A a + v a + Hv a = , (6) 
Q a = cu a - \ (curlt;) a = -\ (curlv) Q , (7) 

Cab = O-ab + D( a U h ) = D( a U fc ) . (8) 
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Notation and some fundamental equations are given in the appendix. The gravito-electric/- 
magnetic Weyl curvature fields are frame-invariant to linear order in v a (cf. the exact boost 
transformation laws (Bll) and (B12) in Jul): 

E a b — E ab , H ab = H ab . (9) 

The gravito-magnetic constraint equation ( A.lOj ), together with (|^) and ([|), shows that 



H ab = = H ab , (10) 
which reflects the absence of tensor modes (gravitational radiation). The div-H ab constraint 



equation (A. 12), together with (j|) and (|l0|), shows that q a is irrotational, and thus so is v a : 

(curlqf = => (curlw) = . (11) 

Thus 

w a = 0, (12) 

by (^), reflecting the absence of vector (rotational) modes. Only scalar modes are admitted, and 
for these 

W a =D a W, (curlWT =0, (13) 
S Qb = D (Q D b) S, (curLS) a6 = 0, (14) 

where W a is a vector orthogonal to u a , S a b is a tracefree tensor orthogonal to u a , and W and S 
are scalars. In particular (to linear order in the second case), 

A a = D Q $ , v a = D a f ^ ) , (15) 



3p, 

where $ is the peculiar gravitational potential, whose existence follows from (0) and the vorticity 



evolution equa tion (A. 4). The velocity potential follows from (|4|) and the div-cr a b constraint 
equation ( |A.9| ). (Note that these potentials have a trivial gauge freedom to add an arbitrary 
homogeneous function.) The velocity potential is exploited in the POTENT method of Bertschinger 



and Dekel |y to "reconstruct" v a from observational data. The shear propagation equation (A. 5) 
then shows that in the present context the gravito-electric Weyl curvature is a purely tidal field 
determined by the peculiar gravitational potential, 

E ab = D (a D b) $ , (16) 

in direct correspondence to the Newtonian situation Q (see also || ) . 

3. PERTURBATION EVOLUTION 

In the (l+3)-covariant approach based on a 4-velocity field u a , because of the expanding Universe 
context, the relativistic equations do not directly provide an analogue of the Newtonian Poisson 
equation for the peculiar gravitational potential. Rather, the equations governing the peculiar 
gravitational potential are integrability conditions, which arise directly from the fact that, with 



Cab = 0, the shear propagation equation (A. 5) is turned into a new constraint equation, i.e. 



Eab = E ab — D{ a A b ) = 0. The time evolution £ a b and the spatial divergence D b £ of this constraint 
equation give us the dynamical equations for the peculiar gravitational potential, which are (in 
units where c = 1 = 8ttG/c 2 ) 0, |l§] 

$ = -±9, (17) 

D a (D 2 $) = iD a/ i-3^D a $-ffD n e ; (18) 

Here D 2 is the covariant Laplace operator for a space of constant curvature K/a 2 , with fC = 0, ±1, 
and a is the background scale factor. Note in particular that ( |l7| ) is not gauge-invariant, but we 
can use it to rewrite the gauge-invariant equation (jlq) as 



D 2 $-3#$ + 3— $ 

a 2 



= D Q [ i/i] , (19) 
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which can be immediately integrated to yield 

D 2 $ - 3ff$ + 3^r$ = i 6(i , (20) 

where we have used the gauge freedom in <I>, i.e., $ — ^ <f> + /3(f), to remove the background part p 
of the matter density p, on the right. This is the (l+3)-covariant relativistic generalization of the 



Newtonian Poisson equation V 2 4> = ^ 5 p. The spatial gradients in (19) ensure gauge-invariance, 
since each term vanishes in the background (note that Sp is gauge-dependent). There are two 
relativistic corrections: one to incorporate background spatial curvature /C, and one to incorporate 
time evolution of the peculiar gravitational potential. The latter is directly induced by the peculiar 
velocity of the matter. [An alternative relativistic Poisson equation is given in the metric-based 
approach by Bardeen [Q, Eq. (4.3).] 

In the linear regime of structure formation, the time scale of variations in $ is of the order of 
the Hubble time H~ l , and on a length scale A, we have D 2 $ ~ $/A 2 . Thus, if A <C H^ 1 , i.e., on 
scales well within the Hubble radius, we have 

\H<Z>\ ~ H 2 \<I>\ < |D 2 $| - M . (21) 

On these scales, where background spatial curvature effects are also negligible, we recover the 
spatial gradient of the Poisson equation. 

3.1. PECULIAR VELOCITY AND PECULIAR ACCELERATION 

Taking the spatial gradient of (fl7|), and using the div-er a fc constraint equation (A. 9), we get 

A a + 2HA a = -\pv a , (22) 



the evolution equation for A a . This is coupled to the Euler momentum equation (A.3) yielding the 
evolution equation for v a as 

v a +Hv a = -A a . (23) 

The source term in (|23|), given by the peculiar acceleration A a = D a $, shows how peculiar velocity 
tends to be generated by the spatial gradient of the peculiar gravitational potential, which is in 
turn generated by matter over- and under-densities through a generalized Poisson-type equation 
derived below. 

It is reasonable to assume that the peculiar velocity is aligned with the peculiar acceleration 
because it is generated by that acceleration , and indeed this assumption is part of the Zel'dovich 
approximation in the weakly non- linear regime p4| . However, in general there will be an orthogonal 
component, i.e., 

v a = -FA a + w a , A a w a = , (24) 
with F a 'growth factor'. Substituting into (p3[), and using (p2[), we get 

w a + (H + \ P F) w a = F — (H — \pF) F - lj A a . (25) 

In order to maintain A a w a = relative to the quasi-Newtonian frame, we find that the 'growth 
factor' F must evolve according to 

F=(H-lpF)F+l+^)w 2 , (26) 

where w 2 = w a w a and A 2 = A a A a . This equation shows, as expected for scalar perturbations, 
that it is impossible to have a purely orthogonal peculiar velocity (F = 0). If we contract j25| ) 
with w a , we find the evolution equation for w. 

w = -2(H+ \pF) w . (27) 

It follows that if w a = at some initial time to, then w a remains zero for all t > t$: irrotational 
gravitational instability does not generate peculiar velocities orthogonal to the peculiar acceleration 
(i.e., to the gradient of the peculiar gravitational potential). If there is an initial orthogonal 
component w a due to random motions, then expansion and matter aggregation serve to decrease 
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it. Furthermore, the rate of change of direction of any initial orthogonal component is always along 
the spatial gradient of writing w a = we a , where e a e a = 1, we find from (p5[) and (p7|) that 



F-{H- \ P F) F - 1 A a . 

We are thus justified in assuming that w a = 0. In this case, 
v a = -FA a , 

where F directly relates v a and A a . Then (^6|) shows that F is a solution of 

F = {H-\pF)F+l , 
while v a follows from (|23|), using (|29|), as 

V Q = . 



V Q exp 



- -H 1 d/ 

i* 1 



(28) 

(29) 
(30) 

(31) 



Here V a is the initial peculiar velocity field, and t is proper time along u a in the background. 
Integration of the background equation ( ^p| ) to find F will thus determine the change with time of 
the peculiar velocity field v a and of the spatial gradient D a <J> = — F~ 1 v a . Then (A. 11) gives the 
corresponding D a /z. 

It is important when solving the Einstein field equations to check that all the field equations 
and their associated consistency conditions are satisfied. This has been verified for the above 



approach in the linearised case in 1 10 



3.2. DENSITY, VELOCITY AND POTENTIAL PERTURBATIONS 

The matter density perturbations are (l+3)-covariantly described by [p 

a 



aB n A a 



(32) 



and are given in the present case by pq | 

A = A + 3a 2 HB a v a , (33) 

where A are the matter density perturbations in the comoving frame. The latter satisfy the 
evolution equation 

(A)"+2JJ(A)"-ipA = , (34) 
which can be solved in a given background. Then the matter density perturbations follow, using 

©,88 

dt~ 



A = A + 3aH(D a V a ) exp 



F 



(35) 



The usual theory of structure formation in the expanding Universe context follows (see ^] , and 
references therein). 

By analogy with the Newtonian relation (2ll [L5l E0| 



v = g , 

p 

we define the position-independent dimcnsionless growth rate / = pF/2H so that 

p 

Note that, by (Jl5|), we can then write 



D„ 



' i e 

<f> H 

/ 3H 



= 



We can rewrite ( |30| ) to become part of the dynamical system 

/' = -(i-in + n A + /)/ + §«, 



(36) 
(37) 

(38) 

(39) 
(40) 
(41) 
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This evolves / together with the background variables and Qa, the latter representing the 
dynamical effects of the cosmological constant. A prime here denotes a derivative with respect to 
the logarithmic dimensionless time var iable N = ln(a/ao) = — ln(l + z), where z is the redshift. 
Note that by the Friedmann equation (|A.15| ), fi/c = 1 — ft — £1a. Equation (|3l| ) now becomes 

^ = V °exp/ g(^U , (42) 



2(1 + z)f 
while ([35]) becomes 

A = A + 3aH(D a V a ) exp 



■ (43) 

We can use numerical solutions for the growth rate / obtained from ( |39| ) - (|4l]) in order to 
evaluate the integrals in the last set of equations, which are the (l+3)-covariant relativistic 
generalizations of well-known relations derived using Newtonian theory. Within the latter context 
a good approximation to / at z — 0, with a cosmological constant, is given by Lahav et al [ p"5[ as 

/owng-a + ifl + ^W. (44) 



70 V 2 



(See Hamilton [ljj for refinements of this approximation.) 

The time-derivative ter m in the generalized Poisson equation (|l^) can be rewritten, using the 



spatial gradient of ( |17j ) and ( A.9), to give 

D a [D 2 $ — 3H 2 (/ — 1 + fi + $] = D a [5 A 4 ] • (45) 
As before, this equation again integrates to give a Poisson-type equation 

D 2 $-3ff 2 C/--fi K )$ = \6n , (46) 

on using the gauge freedom in $ [ $ — > $ + f3(t) ] to remove the background part p of /z on the 
right. The contrast with the Newtonian equation is obvious. 

3.3. SOLUTIONS 

The solution of (Q) may be found using the Green's function in a constant curvature space for 
the inhomogeneous Hclmholtz equation. In the case of a spatially flat background, i.e., when 
flic = 0=1 — — this Green's function is given in ]jj and leads to 

1 [Suit,?') 



*(t,r) = - — / -P^- exp -H(t)y/3m\r- r'\ dV . (47) 
8tt J I r — r \ L J 

With d37j ) and ( ]l5| ) this determines the peculiar velocity field v a . The exponential term represents 
the relativistic correction to the Newtonian formula. Note that the effect of the cosmological 
constant occurs implicitly via the time-dependence of H and /. The way the Newtonian formulae 
work out in practice is shown in the graphs given by Lahav et al [jL5|. 

On scales A = \ f—r'\ well within the Hubble radius, i.e., A <C H^ 1 , (070 shows that 



$ « $Ncwt (1 - V3fH\) . (48) 

For the Einstein-de Sitter case, i.e., when fl^ = = Q\, ( |39| ) has one particular constant 
solution given by / = 1 which corresponds to a growing mode: 

v a = V a a i/2 ( A = A + 3a 5 / 2 iJ (D a V a ) = Ca , (49) 
where C — and we used A oc a. 
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4. CONCLUSIONS 

In this paper we discussed a transparent geometrical framework, based on (l+3)-covariant general 
relativistic methods, to describe peculiar velocities arising from inhomogeneity over cosmologically 
significant distance scales, that also takes into account a cosmological constant. Peculiar velocities 
are defined relative to a preferred frame. In our case, this is defined by the geometrical conditions 
uj a = a ab = imposed on a non-comoving 4- velocity field u a . These relations are always 
true when Bardeen's quasi-Newtonian gauge is used to describe the gravitational dynamics of 
scalar perturbations. In the (l+3)-covariant approach the perturbation equations follow from a 
consistency analysis of the constraint equations. A relativistic analysis provides a theoretical basis 
for determining the peculiar velocity field on cosmologically significant distance scales. Conceivable 
applications of the present framework lie in analyses such as the POTENT project (see B) to 
determine large-scale streaming motions of galaxies within clusters and superclusters. 

A relativistic analysis can also be applied to determine the effect of peculiar velocities of 
matter structures on the temperature anisotropy in the cosmic microwave background radiation, 
requiring a harmonic analysis of the above variables and equations. This was not carried out in 
the present work, but has been treated in a (l+3)-covariant approach by Gebbie et al O, |l2|, [l3| . 
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Appendix A. Basic equations 

The (l+3)-covariant spatial curls are 

(curl V) a = e abc T> b V c , (curl S) ab = e cd{a T> c S b) d , 

where D a is the spatially projected covariant derivative and e abc is the spatially projected 
alternating tensor. The kinematic quantities are O = D a w a (rate of expansion), A a = u a (4- 
acceleration), u a — — i(curlw) Q (vorticity), and a ab — T)^ a u b ) (rate of shear). The dot denotes 
the directional covariant derivative u a V a and the angle brackets denote the projected symmetric 
tracefree part: = [h^ a c h b d — \h cd h ab \S c d- The gravito-electric/-magnetic Weyl curvature 

fields are 

E a b = C ac bdU C U d = E( ab ) , H ab = \e ac dC cd be U e = -ff( a fc) , 

where C abc d is the Weyl curvature tensor. The linearized (l+3)-covariant evolution equations for 
dust (p = ir ab = 0) arc (in any frame, using units c = 1 = 8irG/c 2 , and including scalar, vector 



and tensor modes) |^0[ |l6| 

A + Da<? Q = -9/*, (A.l) 

e-D Q A a = - i6 2 -ip + A, (A.2) 

q a = -4Hq a -pA a , (A.3) 

w a + i(curlA) a = -2Huo a , (A.4) 

o- a h - D( A b ) = -2Ha ab - E ab , (A. 5) 

E ab - (curl H ) ab + ^~D( a <lb) =—31/ E ab - \ p a ab , (A. 6) 

H ab + (curl E) ab = - 3HH ab , (A.7) 
and the linearized (l+3)-covariant constraint equations are 

- Y) a uj a , (A.8) 

= D & a ab - (curl uj) a - f D a 6 + q a , (A.9) 

= (curlCT) ab + B (a uj b) - H ab , (A.10) 

= F> b E ab - ±V a (i + Hq a , (A.ll) 

= V b H ab + ±{cmlq) a ~ pcj a , (A.12) 
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where H = a/a is the background Hubble expansion rate, related to the background value p of fi 

by 



H 



|(p + A)-^5 , 

i £ 

- 2^+72 • 



(A.13) 
(A.14) 



The Friedmann equation (A.13) may be rewritten in terms of the dimensionless density parameters 
as 



1 = Q + fi A + 



(A.15) 
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